Multiomic analysis of HER2-enriched and AR-positive breast carcinoma with apocrine differentiation and an oligometastatic course: a case report.

Breast carcinoma is the most prevalent cancer among women globally. It has variable clinical courses depending on the stage and clinical-biological features. This case report describes a 56-year-old female with invasive breast cancer without estrogen or progesterone receptor expression, with apocrine differentiation, and with no germline variants in the BRCA1 and BRCA2 genes. Throughout the clinical course, the patient exhibited discordant results for HER2 in immunohistochemistry and in situ hybridization. During the second relapse, the disease displayed apocrine microscopic features. The tumor underwent analysis for the androgen receptor, GCDFP-15, RNA-seq, and whole-genome sequencing (WGS) to identify the breast cancer subtype and to characterize the cancer genome. Our bioinformatic analysis revealed 20,323 somatic SNV/Indels, including five mutations in cancer-related genes that are believed to be responsible for the tumor's development. Two of these mutations were found in the PIK3CA and TP53 genes. Furthermore, the tumor tissue exhibited large copy number alterations to the chromosomes, which could impact gene expression through complex mechanisms and contribute to the tumor phenotype. Clustering algorithms applied on RNA-sequencing data categorized this cancer as a HER2+ subtype. The second-line capecitabine chemotherapy treatment is ongoing, and the patient is responding well. Bioinformatic results support the current treatment decision and open the way to further treatments.

iWhale: a computational pipeline based on Docker and SCons for detection and annotation of somatic variants in cancer WES data.

Whole exome sequencing (WES) is a powerful approach for discovering sequence variants in cancer cells but its time effectiveness is limited by the complexity and issues of WES data analysis. Here we present iWhale, a customizable pipeline based on Docker and SCons, reliably detecting somatic variants by three complementary callers (MuTect2, Strelka2 and VarScan2). The results are combined to obtain a single variant call format file for each sample and variants are annotated by integrating a wide range of information extracted from several reference databases, ultimately allowing variant and gene prioritization according to different criteria. iWhale allows users to conduct a complex series of WES analyses with a powerful yet customizable and easy-to-use tool, running on most operating systems (macOs, GNU/Linux and Windows). iWhale code is freely available at https://github.com/alexcoppe/iWhale and the docker image is downloadable from https://hub.docker.com/r/alexcoppe/iwhale.

A high definition picture of somatic mutations in chronic lymphoproliferative disorder of natural killer cells.

The molecular pathogenesis of chronic lymphoproliferative disorder of natural killer (NK) cells (CLPD-NK) is poorly understood. Following the screening of 57 CLPD-NK patients, only five presented STAT3 mutations. WES profiling of 13 cases negative for STAT3/STAT5B mutations uncovered an average of 18 clonal, population rare and deleterious somatic variants per patient. The mutational landscape of CLPD-NK showed that most patients carry a heavy mutational burden, with major and subclonal deleterious mutations co-existing in the leukemic clone. Somatic mutations hit genes wired to cancer proliferation, survival, and migration pathways, in the first place Ras/MAPK, PI3K-AKT, in addition to JAK/STAT (PIK3R1 and PTK2). We confirmed variants with putative driver role of MAP10, MPZL1, RPS6KA1, SETD1B, TAOK2, TMEM127, and TNFRSF1A genes, and of genes linked to viral infections (DDX3X and RSF1) and DNA repair (PAXIP1). A truncating mutation of the epigenetic regulator TET2 and a variant likely abrogating PIK3R1-negative regulatory activity were validated. This study significantly furthered the view of the genes and pathways involved in CLPD-NK, indicated similarities with aggressive diseases of NK cells and detected mutated genes targetable by approved drugs, being a step forward to personalized precision medicine for CLPD-NK patients.

Lack of Viral Load Within Chronic Lymphoproliferative Disorder of Natural Killer Cells: What Is Outside the Leukemic Clone?

Large granular lymphocyte leukemias (LGLL) are sustained by proliferating cytotoxic T cells or NK cells, as happens in Chronic Lymphoproliferative Disorder of Natural Killer cells (CLPD-NK), whose etiology is only partly understood. Different hypotheses have been proposed on the original events triggering NK cell hyperactivation and transformation, including a role of viral agents. In this perspective, we revise the lines of evidence that suggested a pathogenetic role in LGLL of the exposure to retroviruses and that identified Epstein Barr Virus (EBV) in other NK cell leukemias and lymphomas and focus on the contrasting data about the importance of viral agents in CLPD-NK. EBV was detected in aggressive NK leukemias but not in the indolent CLPD-NK, where seroreactivity against HTLV-1 retrovirus envelope BA21 protein antigens has been reported in patients, although lacking clear evidence of HTLV infection. We next present original results of whole exome sequencing data analysis that failed to identify viral sequences in CLPD-NK. We recently demonstrated that proliferating NK cells of patients harbor several somatic lesions likely contributing to sustain NK cell proliferation. Thus, we explore whether "neoantigens" similar to the BA21 antigen could be generated by aberrancies present in the leukemic clone. In light of the literature and new data, we evaluated the intriguing hypothesis that NK cell activation can be caused by retroviral agents located outside the hematopoietic compartment and on the possible mechanisms involved with the prospects of immunotherapy-based approaches to limit the growth of NK cells in CLPD-NK disease.

Somatic mutations in specific and connected subpathways are associated with short neuroblastoma patients' survival and indicate proteins targetable at onset of disease.

Neuroblastoma (NB) is an embryonic malignancy of the sympathetic nervous system with heterogeneous biological, morphological, genetic and clinical characteristics. Although genomic studies revealed the specific biological features of NB pathogenesis useful for new therapeutic approaches, the improvement of high-risk (HR)-NB patients overall survival remains unsatisfactory. To further clarify the biological basis of disease aggressiveness, we used whole-exome sequencing to examine the genomic landscape of HR-NB patients at stage M with short survival (SS) and long survival (LS). Only a few genes, including SMARCA4, SMO, ZNF44 and CHD2, were recurrently and specifically mutated in the SS group, confirming the low recurrence of common mutations in this tumor. A systems biology approach revealed that in the two patient groups, mutations occurred in different pathways. Mutated genes (ARHGEF11, CACNA1G, FGF4, PTPRA, PTK2, ANK3, SMO, NTNG2, VCL and NID2) regulate the MAPK pathway associated with the organization of the extracellular matrix, cell motility through PTK2 signaling and matrix metalloproteinase activity. Moreover, we detected mutations in LAMA2, PTK2, LAMA4, and MMP14 genes, impairing MET signaling, in SFI1 and CHD2 involved in centrosome maturation and chromosome remodeling, in AK7 and SPTLC2, which regulate the metabolism of nucleotides and lipoproteins, and in NALCN, SLC12A1, SLC9A9, which are involved in the transport of small molecules. Notably, connected networks of somatically mutated genes specific for SS patients were identified. The detection of mutated genes present at the onset of disease may help to address an early treatment of HR-NB patients using FDA-approved compounds targeting the deregulated pathways.

Somatic mutations activating Wiskott-Aldrich syndrome protein concomitant with RAS pathway mutations in juvenile myelomonocytic leukemia patients.

The WAS gene product is expressed exclusively in the cytoplasm of hematopoietic cells and constitutional genetic abrogation of WASP leads to Wiskott-Aldrich syndrome (WAS). Moreover, mutational activation of WASP has been associated with X-linked neutropenia. Although studies reported that patients with constitutional WAS mutations affecting functional WASP expression may present juvenile myelomonocytic leukemia (JMML)-like features, confounding differential diagnosis above all in the copresence of mutated RAS, an activating somatic mutation of WASP has not been previously described in JMML patients. In our ongoing studies on JMML genomics, we at first detected a somatic WAS mutation in a major clone found at two consecutive relapses in one of two twins with JMML. Both twins were treated with hematopoietic stem cell transplantation after diagnosis of JMML. The somatic WAS mutation detected here displayed an activating WASP phenotype. Screening of 46 sporadic JMML patients at disease onset for mutations in the same PBD domain of WAS revealed two additional singleton patients carrying minor mutated clones. This is the first study to associate somatically acquired WASP mutations with a hematopoietic malignancy and increases insight in the complexity of the genomic landscape of JMML that shows low recurrent mutations concomitant with general hyperactivation of RAS pathway signaling.

Small RNA Sequencing Uncovers New miRNAs and moRNAs Differentially Expressed in Normal and Primary Myelofibrosis CD34+ Cells.

Myeloproliferative neoplasms (MPN) are chronic myeloid cancers thought to arise at the level of CD34+ hematopoietic stem/progenitor cells. They include essential thrombocythemia (ET), polycythemia vera (PV) and primary myelofibrosis (PMF). All can progress to acute leukemia, but PMF carries the worst prognosis. Increasing evidences indicate that deregulation of microRNAs (miRNAs) might plays an important role in hematologic malignancies, including MPN. To attain deeper knowledge of short RNAs (sRNAs) expression pattern in CD34+ cells and of their possible role in mediating post-transcriptional regulation in PMF, we sequenced with Illumina HiSeq2000 technology CD34+ cells from healthy subjects and PMF patients. We detected the expression of 784 known miRNAs, with a prevalence of miRNA up-regulation in PMF samples, and discovered 34 new miRNAs and 99 new miRNA-offset RNAs (moRNAs), in CD34+ cells. Thirty-seven small RNAs were differentially expressed in PMF patients compared with healthy subjects, according to microRNA sequencing data. Five miRNAs (miR-10b-5p, miR-19b-3p, miR-29a-3p, miR-379-5p, and miR-543) were deregulated also in PMF granulocytes. Moreover, 3'-moR-128-2 resulted consistently downregulated in PMF according to RNA-seq and qRT-PCR data both in CD34+ cells and granulocytes. Target predictions of these validated small RNAs de-regulated in PMF and functional enrichment analyses highlighted many interesting pathways involved in tumor development and progression, such as signaling by FGFR and DAP12 and Oncogene Induced Senescence. As a whole, data obtained in this study deepened the knowledge of miRNAs and moRNAs altered expression in PMF CD34+ cells and allowed to identify and validate a specific small RNA profile that distinguishes PMF granulocytes from those of normal subjects. We thus provided new information regarding the possible role of miRNAs and, specifically, of new moRNAs in this disease.

mRNA-Seq and microarray development for the Grooved Carpet shell clam, Ruditapes decussatus: a functional approach to unravel host-parasite interaction.

BACKGROUND: The Grooved Carpet shell clam Ruditapes decussatus is the autochthonous European clam and the most appreciated from a gastronomic and economic point of view. The production is in decline due to several factors such as Perkinsiosis and habitat invasion and competition by the introduced exotic species, the manila clam Ruditapes philippinarum. After we sequenced R. decussatus transcriptome we have designed an oligo microarray capable of contributing to provide some clues on molecular response of the clam to Perkinsiosis. RESULTS: A database consisting of 41,119 unique transcripts was constructed, of which 12,479 (30.3%) were annotated by similarity. An oligo-DNA microarray platform was then designed and applied to profile gene expression in R. decussatus heavily infected by Perkinsus olseni. Functional annotation of differentially expressed genes between those two conditionswas performed by gene set enrichment analysis. As expected, microarrays unveil genes related with stress/infectious agents such as hydrolases, proteases and others. The extensive role of innate immune system was also analyzed and effect of parasitosis upon expression of important molecules such as lectins reviewed. CONCLUSIONS: This study represents a first attempt to characterize Ruditapes decussatus transcriptome, an important marine resource for the European aquaculture. The trancriptome sequencing and consequent annotation will increase the available tools and resources for this specie, introducing the possibility of high throughput experiments such as microarrays analysis. In this specific case microarray approach was used to unveil some important aspects of host-parasite interaction between the Carpet shell clam and Perkinsus, two non-model species, highlighting some genes associated with this interaction. Ample information was obtained to identify biological processes significantly enriched among differentially expressed genes in Perkinsus infected versus non-infected gills. An overview on the genes related with the immune system on R. decussatus transcriptome is also reported.

Transcriptome sequencing and de novo annotation of the critically endangered Adriatic sturgeon.

BACKGROUND: Sturgeons are a group of Condrostean fish with very high evolutionary, economical and conservation interest. The eggs of these living fossils represent one of the most high prized foods of animal origin. The intense fishing pressure on wild stocks to harvest caviar has caused in the last decades a dramatic decline of their distribution and abundance leading the International Union for Conservation of Nature to list them as the more endangered group of species. As a direct consequence, world-wide efforts have been made to develop sturgeon aquaculture programmes for caviar production. In this context, the characterization of the genes involved in sex determination could provide relevant information for the selective farming of the more profitable females. RESULTS: The 454 sequencing of two cDNA libraries from the gonads and brain of one male and one female full-sib A. naccarii, yielded 182,066 and 167,776 reads respectively, which, after strict quality control, were iterative assembled into more than 55,000 high quality ESTs. The average per-base coverage reached by assembling the two libraries was 4X. The multi-step annotation process resulted in 16% successfully annotated sequences with GO terms. We screened the transcriptome for 32 sex-related genes and highlighted 7 genes that are potentially specifically expressed, 5 in male and 2 in females, at the first life stage at which sex is histologically identifiable. In addition we identified 21,791 putative EST-linked SNPs and 5,295 SSRs. CONCLUSIONS: This study represents the first large massive release of sturgeon transcriptome information that we organized into the public database AnaccariiBase, which is freely available at http://compgen.bio.unipd.it/anaccariibase/. This transcriptomic data represents an important source of information for further studies on sturgeon species. The hundreds of putative EST-linked molecular makers discovered in this study will be invaluable for sturgeon reintroduction and breeding programs.

Genome evolution in the cold: Antarctic icefish muscle transcriptome reveals selective duplications increasing mitochondrial function.

Antarctic notothenioids radiated over millions of years in subzero waters, evolving peculiar features, such as antifreeze glycoproteins and absence of heat shock response. Icefish, family Channichthyidae, also lack oxygen-binding proteins and display extreme modifications, including high mitochondrial densities in aerobic tissues. A genomic expansion accompanying the evolution of these fish was reported, but paucity of genomic information limits the understanding of notothenioid cold adaptation. We reconstructed and annotated the first skeletal muscle transcriptome of the icefish Chionodraco hamatus providing a new resource for icefish genomics (http://compgen.bio.unipd.it/chamatusbase/, last accessed December 12, 2012). We exploited deep sequencing of this energy-dependent tissue to test the hypothesis of selective duplication of genes involved in mitochondrial function. We developed a bioinformatic approach to univocally assign C. hamatus transcripts to orthology groups extracted from phylogenetic trees of five model species. Chionodraco hamatus duplicates were recorded for each orthology group allowing the identification of duplicated genes specific to the icefish lineage. Significantly more duplicates were found in the icefish when transcriptome data were compared with whole-genome data of model species. Indeed, duplicated genes were significantly enriched in proteins with mitochondrial localization, involved in mitochondrial function and biogenesis. In cold conditions and without oxygen-carrying proteins, energy production is challenging. The combination of high mitochondrial densities and the maintenance of duplicated genes involved in mitochondrial biogenesis and aerobic respiration might confer a selective advantage by improving oxygen diffusion and energy supply to aerobic tissues. Our results provide new insights into the genomic basis of icefish cold adaptation.

Sequencing and characterization of striped venus transcriptome expand resources for clam fishery genetics.

BACKGROUND: The striped venus Chamelea gallina clam fishery is among the oldest and the largest in the Mediterranean Sea, particularly in the inshore waters of northern Adriatic Sea. The high fishing pressure has lead to a strong stock abundance decline, enhanced by several irregular mortality events. The nearly complete lack of molecular characterization limits the available genetic resources for C. gallina. We achieved the first transcriptome of this species with the aim of identifying an informative set of expressed genes, potential markers to assess genetic structure of natural populations and molecular resources for pathogenic contamination detection. METHODOLOGY/PRINCIPAL FINDINGS: The 454-pyrosequencing of a normalized cDNA library of a pool C. gallina adult individuals yielded 298,494 raw reads. Different steps of reads assembly and filtering produced 36,422 contigs of high quality, one half of which (18,196) were annotated by similarity. A total of 111 microsatellites and 20,377 putative SNPs were identified. A panel of 13 polymorphic transcript-linked microsatellites was developed and their variability assessed in 12 individuals. Remarkably, a scan to search for contamination sequences of infectious origin indicated the presence of several Vibrionales species reported to be among the most frequent clam pathogen's species. Results reported in this study were included in a dedicated database available at http://compgen.bio.unipd.it/chameleabase. CONCLUSIONS/SIGNIFICANCE: This study represents the first attempt to sequence and de novo annotate the transcriptome of the clam C. gallina. The availability of this transcriptome opens new perspectives in the study of biochemical and physiological role of gene products and their responses to large and small-scale environmental stress in C. gallina, with high throughput experiments such as custom microarray or targeted re-sequencing. Molecular markers, such as the already optimized EST-linked microsatellites and the discovered SNPs will be useful to estimate effects of demographic processes and to detect minute levels of population structuring.

Surviving in a toxic world: transcriptomics and gene expression profiling in response to environmental pollution in the critically endangered European eel.

BACKGROUND: Genomic and transcriptomic approaches have the potential for unveiling the genome-wide response to environmental perturbations. The abundance of the catadromous European eel (Anguilla anguilla) stock has been declining since the 1980s probably due to a combination of anthropogenic and climatic factors. In this paper, we explore the transcriptomic dynamics between individuals from high (river Tiber, Italy) and low pollution (lake Bolsena, Italy) environments, which were measured for 36 PCBs, several organochlorine pesticides and brominated flame retardants and nine metals. RESULTS: To this end, we first (i) updated the European eel transcriptome using deep sequencing data with a total of 640,040 reads assembled into 44,896 contigs (Eeelbase release 2.0), and (ii) developed a transcriptomic platform for global gene expression profiling in the critically endangered European eel of about 15,000 annotated contigs, which was applied to detect differentially expressed genes between polluted sites. Several detoxification genes related to metabolism of pollutants were upregulated in the highly polluted site, including genes that take part in phase I of the xenobiotic metabolism (CYP3A), phase II (glutathione-S-transferase) and oxidative stress (glutathione peroxidase). In addition, key genes in the mitochondrial respiratory chain and oxidative phosphorylation were down-regulated at the Tiber site relative to the Bolsena site. CONCLUSIONS: Together with the induced high expression of detoxification genes, the suggested lowered expression of genes supposedly involved in metabolism suggests that pollution may also be associated with decreased respiratory and energy production.

MAGIA²: from miRNA and genes expression data integrative analysis to microRNA-transcription factor mixed regulatory circuits (2012 update).

MAGIA(2) (http://gencomp.bio.unipd.it/magia2) is an update, extension and evolution of the MAGIA web tool. It is dedicated to the integrated analysis of in silico target prediction, microRNA (miRNA) and gene expression data for the reconstruction of post-transcriptional regulatory networks. miRNAs are fundamental post-transcriptional regulators of several key biological and pathological processes. As miRNAs act prevalently through target degradation, their expression profiles are expected to be inversely correlated to those of the target genes. Low specificity of target prediction algorithms makes integration approaches an interesting solution for target prediction refinement. MAGIA(2) performs this integrative approach supporting different association measures, multiple organisms and almost all target predictions algorithms. Nevertheless, miRNAs activity should be viewed as part of a more complex scenario where regulatory elements and their interactors generate a highly connected network and where gene expression profiles are the result of different levels of regulation. The updated MAGIA(2) tries to dissect this complexity by reconstructing mixed regulatory circuits involving either miRNA or transcription factor (TF) as regulators. Two types of circuits are identified: (i) a TF that regulates both a miRNA and its target and (ii) a miRNA that regulates both a TF and its target.

Impact of host genes and strand selection on miRNA and miRNA* expression.

Dysregulation of miRNAs expression plays a critical role in the pathogenesis of genetic, multifactorial disorders and in human cancers. We exploited sequence, genomic and expression information to investigate two main aspects of post-transcriptional regulation in miRNA biogenesis, namely strand selection regulation and expression relationships between intragenic miRNAs and host genes. We considered miRNAs expression profiles, measured in five sizeable microarray datasets, including samples from different normal cell types and tissues, as well as different tumours and disease states. First, the study of expression profiles of "sister" miRNA pairs (miRNA/miRNA*, 5' and 3' strands of the same hairpin precursor) showed that the strand selection is highly regulated since it shows tissue-/cell-/condition-specific modulation. We used information about the direction and the strength of the strand selection bias to perform an unsupervised cluster analysis for the sample classification evidencing that is able to distinguish among different tissues, and sometimes between normal and malignant cells. Then, considering a minimum expression threshold, in few miRNA pairs only one mature miRNA is always present in all considered cell types, whereas the majority of pairs were concurrently expressed in some cell types and alternatively in others. In a significant fraction of concurrently expressed pairs, the major and the minor forms found at comparable levels may contribute to post-transcriptional gene silencing, possibly in a coordinate way. In the second part of the study, the behaved tendency to co-expression of intragenic miRNAs and their "host" mRNA genes was confuted by expression profiles examination, suggesting that the expression profile of a given host gene can hardly be a good estimator of co-transcribed miRNA(s) for post-transcriptional regulatory networks inference. Our results point out the regulatory importance of post-transcriptional phases of miRNAs biogenesis, reinforcing the role of such layer of miRNA biogenesis in miRNA-based regulation of cell activities.

Transcriptome sequencing and microarray development for the Manila clam, Ruditapes philippinarum: genomic tools for environmental monitoring.

BACKGROUND: The Manila clam, Ruditapes philippinarum, is one of the major aquaculture species in the world and a potential sentinel organism for monitoring the status of marine ecosystems. However, genomic resources for R. philippinarum are still extremely limited. Global analysis of gene expression profiles is increasingly used to evaluate the biological effects of various environmental stressors on aquatic animals under either artificial conditions or in the wild. Here, we report on the development of a transcriptomic platform for global gene expression profiling in the Manila clam. RESULTS: A normalized cDNA library representing a mixture of adult tissues was sequenced using a ultra high-throughput sequencing technology (Roche 454). A database consisting of 32,606 unique transcripts was constructed, 9,747 (30%) of which could be annotated by similarity. An oligo-DNA microarray platform was designed and applied to profile gene expression of digestive gland and gills. Functional annotation of differentially expressed genes between different tissues was performed by enrichment analysis. Expression of Natural Antisense Transcripts (NAT) analysis was also performed and bi-directional transcription appears a common phenomenon in the R. philippinarum transcriptome. A preliminary study on clam samples collected in a highly polluted area of the Venice Lagoon demonstrated the applicability of genomic tools to environmental monitoring. CONCLUSIONS: The transcriptomic platform developed for the Manila clam confirmed the high level of reproducibility of current microarray technology. Next-generation sequencing provided a good representation of the clam transcriptome. Despite the known limitations in transcript annotation and sequence coverage for non model species, sufficient information was obtained to identify a large set of genes potentially involved in cellular response to environmental stress.

TRAM (Transcriptome Mapper): database-driven creation and analysis of transcriptome maps from multiple sources.

BACKGROUND: Several tools have been developed to perform global gene expression profile data analysis, to search for specific chromosomal regions whose features meet defined criteria as well as to study neighbouring gene expression. However, most of these tools are tailored for a specific use in a particular context (e.g. they are species-specific, or limited to a particular data format) and they typically accept only gene lists as input. RESULTS: TRAM (Transcriptome Mapper) is a new general tool that allows the simple generation and analysis of quantitative transcriptome maps, starting from any source listing gene expression values for a given gene set (e.g. expression microarrays), implemented as a relational database. It includes a parser able to assign univocal and updated gene symbols to gene identifiers from different data sources. Moreover, TRAM is able to perform intra-sample and inter-sample data normalization, including an original variant of quantile normalization (scaled quantile), useful to normalize data from platforms with highly different numbers of investigated genes. When in 'Map' mode, the software generates a quantitative representation of the transcriptome of a sample (or of a pool of samples) and identifies if segments of defined lengths are over/under-expressed compared to the desired threshold. When in 'Cluster' mode, the software searches for a set of over/under-expressed consecutive genes. Statistical significance for all results is calculated with respect to genes localized on the same chromosome or to all genome genes. Transcriptome maps, showing differential expression between two sample groups, relative to two different biological conditions, may be easily generated. We present the results of a biological model test, based on a meta-analysis comparison between a sample pool of human CD34+ hematopoietic progenitor cells and a sample pool of megakaryocytic cells. Biologically relevant chromosomal segments and gene clusters with differential expression during the differentiation toward megakaryocyte were identified. CONCLUSIONS: TRAM is designed to create, and statistically analyze, quantitative transcriptome maps, based on gene expression data from multiple sources. The release includes FileMaker Pro database management runtime application and it is freely available at http://apollo11.isto.unibo.it/software/, along with preconfigured implementations for mapping of human, mouse and zebrafish transcriptomes.